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ABSTRACT 

A numerical study of incompressible turbulent separated flows is carried out by 
using two-equation turbulence models of the K-e type. On the basis of realizabil- 
ity analysis, a new formulation of the eddy-viscosity is proposed which ensures the 
positiveness of turbulent normal stresses — a realizability condition that most ex- 
isting two-equation turbulence models are unable to satisfy. The present model is 
applied to calculate two backward-facing step flows. Calculations with the stan- 
dard K-e model and a recently developed RNG-based K-e model are also made for 
comparison. The calculations are performed with a finite-volume method. A second- 
order accurate differencing scheme and sufficiently fine grids are used to ensure the 
numerical accuracy of solutions. The calculated results are compared with the ex- 
perimental data for both mean and turbulent quantities. The comparison shows 
that the present model performs quite well for separated flows. 


1. INTRODUCTION 


Turbulent separated flows occur in a number of engineering applications. Because 
of their great practical importance, there is a strong demand for calculation methods 
to predict such flows. Turbulent flow over a backward- facing step is one of the most 
extensively used benchmark cases in the study of turbulence models for separated 
flows. It involves severe adverse pressure gradient, streamline curvature, coexistence 
of both strong and weak shear layers as well as significant extra strain rates in more 
than one direction, thereby constituting a severe test for turbulence models. If a 
turbulence model can correctly simulate this flow, it will be likely to be successful 
with other complicated flows. 

The relevant experimented studies on backward-facing step flows are reported in 
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Bradshaw and Wong (1972), Driver and SeegmiUer (1985), Driver el al. (1987), 
Durst and Schmitt (1985), Eaton and Johnston (1980), Kim et al. (1978, 1980), 
Stevenson et al. (1984) and Westphal et al. (1981). Among them, the case of Kim 
et al. (1978) with a larger expansion was a test case (0421) for the 1980-81 Stanford 
Conference on Complex Turbulent Flows (Kline et al., 1981), which has extensively 
been used to validate numerical calculations. However, this case has no turbulent 
data in the recirculation zone. The case of Driver and Seegmiller (1985) with a 
smaller expansion provides detailed data, including the wall friction coefficient and 
the turbulent quantities up to triple correlations. 

The recent calculations with turbulence modeling can be found in Avva et al. 
(1990), Celenligil and Mellor (1985), Obi et al. (1989), So and Lai (1988), Speziale 
and Ngo (1988), Speziale and Thangam (1992) and Thangam and Hur (1991). The 
calculations of Celenligil and Mellor, Obi et al., and So and Lai were carried out 
with second-order closures, and the others with the standard K-e model and its 
variants. These calculations show that the K-e model largely underpredicts the 
reattachment point which is a sensitive parameter to assess the overall performance 
of turbulence models. No definitive conclusion can be drawn with the second-order 
closures, because Celenligil and Mellor obtained an overprediction, while Obi et 
al. and So and Lai obtained an underprediction of the reattachment point. The 
overall improvement achieved with these second-order closures is not strong enough 
to establish their convincing superiority over the K-e model in calculating separated 
flows. 

In the standard K-e model, all the model coefficients are constant which are de- 
termined from a set of experiments for simple turbulent flows. Numerical experience 
over the last two decades has shown that this set of constants have a broad appli- 
cability, but they should not be expected to be universal. Rodi (1972) found that 
the K-e model’s ability to predict weak shear flows can be significantly improved 
by using as a function of the average ratio of P/e (P is the production of the 
turbulent kinetic energy) instead of a constant. Leschziner and Rodi (1981) pro- 
posed a function for which takes into account th e effect o f streamline curvature 
and obtained improved results in the calculation of annular and twin parallel jets. 
Recently, Yakhot and co-workers have developed a version of the K-e model using 
Renormalization Group (RNG) method. This model is of the same form as the stan- 
dard -K-e model, but all the model coefficients assume different values. In the latest 
version of the RNG based K-e model (Speziale and Thangam, 1992), the coefficient 
Ci related to the production of dissipation term is set to a function of 77 , where 77 
is the time scale ratio of the turbulence to the mean flow field. The reattachment 
point predicted by this model is within 5% of the experimental value for the case of 
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Kim et al. (1978). 

In this study, the realizability principle (Schumann, 1977 and Lumley, 1978) is 
applied to analyze the K-e model. The analysis results in a new formulation of 
which is a function of time scale ratio of the turbulence to the mean strain rate. 
The new C7 M will ensure the positivity of each component of the turbulent kinetic 
energy - realizability that most existing eddy-viscosity models do not satisfy. The 
model validation is made on the basis of applications to the two backward-facing 
step flows experimently studied by Driver and Seegmiller (1985) and Kim et al. 
(1978). Calculations are carried out with a conservative finite- volume method, and 
a second-order accurate and bounded differencing scheme together with sufficiently 
fine grids is used to ensure the solution both grid-independent and free from nu- 
merical diffusion. The calculated results are compared in detail with experimental 
data as well as with those obtained using the standard K-e model and the RNG K-e 
model. 


2. MATHEMATICAL FORMULATION 
2.1 Governing Equations 

For incompressible steady flows, the non-dimensional governing equations formu- 
lated wit hin the framework of the K-e model may be written as 


u jtj = 0 

( 1 ) 

(UjUi - + u i u i)j = -P.i 

( 2 ) 


( 3 ) 

M‘-(Te + 7. )€ ^ = c 'ic p - C 4 

( 4 ) 

- UiUj = -~KSij + u t (Uij + Uj t i) 

( 5 ) 

CM j 

o* 

II 

£ 

( 6 ) 
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P = -UiUjUij (7) 

where non-dimensionalization is made by using the reference length L re f and the 
reference velocity U re f • Accordingly, the flow Reynolds number is defined by 

Re = hl£?L ( 8 ) 

V 

In the standard K-e model (Launder and Spalding, 1974), the model coefficients 
Cl, C 2 , <Jk and <Tf assume the following constant values: 

= 0.09, Ci = 1.44, C 2 = 1.92, a K = 1, <r 4 = 1.3 (9) 

and in the RNG K-e model (Speziale and Thangam, 1992), they are: 


= 0.085, Ci - 1.42 - 


7?(1 — 77/ 4.38) 
1 + 0.0157/ 3 


, C 2 = 1.68, <t k = C7-* = 0.7179 (10) 


where 


and 


i) = SK/e, S = (2SijSij) 1/2 




(ID 

( 12 ) 


2.2 Realizability 

Realizability (Schumann, 1977, Lumley, 1978) which requires the non- negativity 
of turbulent normal stresses is a basic physical and mathematical principle that the 
solution of any turbulence model equation should obey. It also represents a minimal 
requirement to prevent a turbulence model from producing unphysical results. In the 
following, we will apply this principle to derive constraint on the model coefficients. 
Consider a deformation rate tensor of the form 

( Ui,i 0 0 \ 

0 C 2 , 2 0 (13) 

V 0 00/ 

The continuity equation ( 1) gives 

U 2>2 = —U\ t \ (14) 
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and from Eq. ( 5 ), the normal stress UjUi can be written as 


U\U\ 

T 


2 

3 


-c*n 


Note that in case of Eqs. ( 13 ) and ( 14 ), 77 can be written as 


( 15 ) 


77 = 


2Ui,iK 

e 


( 16 ) 


Physically, u^u{ will decrease with an increase in the mean strain rate Ui t u but 
cannot be driven to negative values. Therefore, realizability conditions for 

are: 


U\U\ 

K 

> 0, if 0 < 77 < oo 

( 17 ) 

UlUi 

i 

0 

4 

■3 

1 
8 

( 18 ) 

K 



,UlUi v 

{ K ),v 

— i 0, if 77 -+ 00 

( 19 ) 

These conditions can be satisfied 

by specifying as: 




(20) 


M A + 77 


where A is a positive constant. 

Similar analysis on also leads to Eq. ( 20 ). It should be mentioned that Eq. 
( 20) also holds in the case of three-dimensional pure strain rates 


u ltl 0 0 \ 

0 Uz,2 o 

0 0 Uz t z ) 


( 21 ) 


and that any deformation rate tensor can be written in the form of ( 21) in the 

principal axes of deformation rate tensor. 

The use of Eq.( 20 ) while keeping the other model coefficients the same as those 
in the standard K-e model constitutes the present realizable isotropic K-e model. 
The value of the extra model constant A is taken as 


A = 5.5 


( 22 ) 


which has been found to work well for both the test cases considered in this study. 
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3. NUMERICAL SOLUTION 


In two dimensions, the transport equations ( 1) to ( 4) can be written in the 
following general form 

[U\4> - (-^ + ^)4>,i],i + [U2<f> - (-^ + —)<!>, 2], 2 = S4, (23) 

where <f> stands for U U U 2 ,K or e. For the momentum equations, the source term 
Sj, includes the cross- derivative diffusion terms. 

The numerical method used to solve the system of equations ( 23) is a finite- 
volume procedure. It uses a non-staggered grid with all the dependent variables 
being stored at the same geometric center of each control volume. The momentum 
interpolation procedure of Rhie and Chow (1983) is used to avoid spurious oscilla- 
tions usually associated with the non-staggered grid, and the pressure-velocity cou- 
pling is handled with the SIMPLEC algorithm (Van Doormal and Raithby, 1984). 
To ensure both accuracy and stability of numerical solution, the convection terms are 
approximated by a second-order and bounded differencing scheme (Zhu, 1991a), and 
all the other terms by the conventional central differencing scheme. The strongly 
implicit procedure of Stone (1968) is used to solve Jhe system of algebraic equa- 
tions. The iterative solution process is considered converged when the maximum 
normalised residue of all the dependent variables is less than 10 . The details of 

the present numerical procedure are given in Zhu (1991b). 


4. RESULTS AND DISCUSSION 

The present model together with the standard K-e model and the RNG K-e model 
are applied to the two back war d-facing step flows experimentally studied by Kim, 
Kline and Johnston (1978) and Driver and Seegnnller (1985), from here on referred 
to as KKJ- and DS-cases, respectively. Fig.l shows the flow configuration and the 
Cartesian co-ordinate system used. Table 1 gives the flow parameters for both cases; 
here the experimental reference free-stream velocities and step heights are taken as 
the reference quantities for non-dimensionalization. 

Table 1. Flow parameters 


case 

Re 

6 

L, 

U 

H. 

H d 

u ref 

KKJ 

44737 

0.6 

10 

40 

1 

2 

1 

DS 

37423 

1.5 

10 

40 

1 

8 

1 
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Boundaries of the flows axe inlet, outlet and solid wall. At the inlet, the ex- 
perimental data are available for the streamwise mean velocity U and the turbulent 
normal stresses uu and vv. K is calculated from these uu and vv with the assumption 
that 


and e by 


ww = —(uu + vv) 

2 V 


(24) 


c 3 J a k 3 > 2 


L = min(0.4lAy, 0.0855) 


(25) 


where A y is the distance from the wall and 8 is the boundary-layer thickness given 
in Table 1. At the outlet, the streamwise derivatives of the flow variables are set 
to zero. Influences of both inlet and outlet conditions on the solution are examined 
by ch an ging the locations L, and £ e , and it has been found that in both cases, 
the distances given in Table 1 are already sufficiently fax away from the region of 
interest. The standard wall function approach (Launder and Spalding, 1974) is used 
to bridge the viscous sublayer near the wall. 

Grid dependence of solutions is examined by using two sets of non-uniform numer- 
ical grids which contain 110x52 (coarse) and 199x91 (fine) points for the KKJ-case 
and 106x56 (coarse) and 201x109 (fine) points for the DS-case. Fig.2(a) shows 
the friction coefficient Cf at the bottom wall, calculated with the present model 
on the two grids in the KKJ-case. It can be seen that the grid refinement from 
110x52 to 199x91 points does produce a noticeable difference. The same also holds 
true for the other two models. This indicates that the solutions obtained on the 
coarse grids have not yet been sufficiently close to the grid-independent stage. Re- 
cently, Thangam and Hur (1991) have conducted a highly-resolved calculation in the 
KKJ-case. They have found that quadrupling a 166x73 grid leads to only a min- 
imal improvement. Therefore, the results with the fine grids can be considered as 
grid-independent. In the DS-case, the fine grid computations required 681/766/800 
iterations and took approximately 8. 3/9. 3/9. 8 minutes of CPU time for the stan- 
dard/RNG/present model on the Cray YMP computer. Only find grid results will 
be presented in the following. 

In Fig. 2(b) the calculated friction coefficients with the three models are compared 
with the experimental data in the DS-case. No such experimental data are available 
in the KKJ-case. It can be seen from Fig.2(b) that all the three models largely un- 
derpredict the negative peak of Cf, pointing to limited accuracy of the wall function 
approach in the recirculation region. In the recovery region and downwards, the 
standard K-e model agrees well with the experimental data, while both the RNG 
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and the present models basically give the same results which are somewhat under- 
predicted. For lack of good near-wall turbulence models for separated flows, it is 
difficult to judge the performance of the models with Cf that is very sensitive to the 
near-wall turbulence modeling. 

Table 2 compares the computed and measured reattachment points. They are 
determined in the calculation from the point where (7/ goes to zero. The reattach- 
ment point is a critical parameter which has often been used to assess the overall 
performance of turbulence models. Table 2 clearly shows that the results of both 
the present and the RNG models are much better than those of the standard model. 


Table 2. Comparison of reattachment points 


case 

experiment 

standard 

RNG 

present 

KKJ 

7 ±0.5 

6.35 

7.47 

7.34 

DS 

6.1 

4.99 

6.01 

5.77 


Figs. 3(a) and 3(b) show the comparison of computed and measured static pressure 
coefficient C p along the bottom waill. In both cases, the standard K-e model is seen 
to predict premature pressure rises, which is consistent with its underprediction 
of the reattachment lengths, while both the present and the RNG models capture 
these pressure rises quite well. The results of both the present and the RNG models 
are very similar, and only at the lower end of steep gradients can some noticeable 
difference be seen. 

The streamwise mean velocity U profiles are shown in Figs.4(a) and 4(b) at four 
different downstream locations. Here again, the present and the RNG models yield 
essentially the same results. They predict reverse flows better than the standard 
K-e model, but result in somewhat slower recovery in regions near the reattachment 
point. Interestingly enough, such a slower recovery has also been found in the RSM 
prediction by Obi et al. (1989). Further downstream, say at x=20 in Fig.4(b), the 
results of the three models nearly coincide with each other. 

In the KKJ-experiment, a high degree of flow unsteadiness was present, causing 
the reattachment point to swing constantly within a range of one step height. As 
a result, no experimental data for turbulent quantities were available in the recir- 
culation region. Conversely, the DS-experiment showed a lower unsteadiness of the 
flow and a smaller uncertainty of the reattachment location. Detailed turbulent 
data were provided in the whole -region of interest. Therefore, the comparison of 
turbulent quantities are restricted only to the DS-case. Figs. 5 and 6 show the com- 
parison of predicted and measured turbulent stresses at four x-lo cations, two before 
and two after the reattachment point. It is seen from Fig.5 that the standard K-e 
model overpredicts the turbulent shear stress all along the flow region, while the 
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present and the RNG models give a better agreement with the experimental data. 
The results of the present and the RNG models are virtually the same except in 
the near-step region (x=2) where the RNG model gives a large underprediction. 
For the turbulent normal stresses in Fig.6, the RNG profiles differ from the present 
profiles. The RNG model largely underpredicts the turbulent normal stresses in the 
recirculation region (x=2 and 5). The present model produces the best results of 
all. These different results of the models may be traced to the different levels of 
the turbulent eddy-viscosity they predict. Fig. 7 shows the turbulent eddy-viscosity 
profiles of the three models in the DS-case. The present and the RNG models con- 
siderably reduce the value of v t , but this reduction is more than enough for the RNG 
model in the near-step region (x=2), resulting in the large underpredictions of the 
turbulent stresses there. 


5. CONCLUSIONS 

A new version of the K-c model has been developed in which the model coefficient 
Op is related to the time scale ratio of the turbulence to the mean strain rate 
through the realizability analysis. The new model ensures the positivity of individual 
turbulent normal stresses, while the standard K-e model, like many others, can 
only ensure the positiveness of the turbulent kinetic energy — sum of the turbulent 
normal stresses. The present model has been compared with the standard K-e model 
and the recently proposed RNG K-e model as well as with the experiments in the 
calculations of the two backward- facing step flows. The comparison shows that 
the present model effectively reduces the turbulent eddy-viscosity level, resulting in 
significant improvement over the standard K-e model. The RNG model generally 
gives very s imil ar predictions to the present model, but overly reduces the turbulent 
eddy- viscosity level in the recirculation region near the step. It should be noted that 
the set of model constants in the standard K-e model have a broad generality and 
have stood the test of time. The present model differs from the standard K-e model 
only in one model coefficient, while all the model coefficients in the RNG model are 
different from the standard values. Therefore, the present model could be expected 
to be more general than the RNG model. 
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Figure 1. Backward-facing step geometry 
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Figure 3. Static pressure coefficient C p along the bottom wall 
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Figure 5. Turbulent shear stress profiles (key to symbols as in figure 3) 
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Figure 6. Turbulent normal stress profiles (key to symbols as in figure 3) 
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Figure 7. Turbulent eddy-viscosity profiles (key to symbols as in figure 3) 
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